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We present results of inelastic neutron scattering study on two bond disordered quasi two- 
dimensional quantum magnets (C4Hi2N2)Cu2(Cli-ccBrcc)6 with x=0.035 and 0.075. We observe 
the increase of spin gap, reduction of magnon bandwidth and a decrease of magnon lifetimes com- 
pared to x=0 sample. Additional magnon damping is observed at higher energies away from zone 
center which is found to follow the density of single particle states. 



I. INTRODUCTION 

Disordered quantum magnets are expected to ex- 
hibit qualitatively different behaviour from the pure sys- 
tems. Apart from straightforward chemical substitu- 
tion of the magnetic ions for non-magnetic counterparts 
a more recently implemented idea focuses on creating 
disorder on peripheral sites involved in superexchange 
interact ions. ^2^5] predictions and scarce exper- 

imental work dealing with stability of spin gaji^, lifetime 
of excitations,^ ^ appearance of in-gap localized states^ 
and novel quantum phases^ ^ in the presence of disorder 
necessitate further experiments. In current paper we dis- 
cuss the impact of random bond disorder on spin dynam- 
ics in a quasi-two-dimensional gapped quantum magnet 
as seen by inelastic neutron scattering experiments. Our 
choice of parent compound is a well-characterized Heisen- 
berg quantum magnet piperazinium hexachlorodicuprate 
(PHCC).^^It features a layered structure with com- 
plex spin network of S=l/2 Cu^+ ions bridged by Cu-Cl- 
Cl-Cu superexchange pathways, see Fig{l] The ground 
state is a nonmagnetic spin liquid and the low energy ex- 
citation spectrum consists of a sharp magnon mode with 
a gap A = 1 meV and a bandwidth of about 1.7 meV. 
Considerable attention has been devoted to magnon life- 
times in PHCC due to observed single magnon decay 
into two magnons as the two energy scales meet."^^^^ 
Our focus is on determining the magnon spectrum, it's 
temperature dependence and magnon lifetimes in bond 
disordered PHCC derivates. We find that although no 
drastic modification of spectrum occurs bond disorder 
increases the spin gap, reduces the magnon bandwidth 
and lifetime. A random phase approximation treatment 
of an effective dimer model can account for temperature 
induced blueshift of the gap and drop of intensity in this 
complex material. The most notable effect of bond disor- 
der is increased magnon lifetimes which can be correlated 
with the single-magnon density of states function. 



II. SAMPLES AND EXPERIMENTAL 
METHODS 

PHCC crystallizes in the triclinic PI group with lat- 
tice parameters a = 7.9691 (5) A, b = 7.0348(5) A, 



c = 6.0836(4)A and angles a = 111.083(3)°, /3 = 
99.947(3)°, 7 = 81.287(4)°. Bond disordered PHCC with 
3.5% and 7.5% nominal Br content (referred to as PHCX) 
were prepared by adjusting the HChHBr solvent ratio 
in the starting solutions. Large fully deuterated sin- 
gle crystals were grown by thermal gradient method.!^ 
Samples used in inelastic neutron scattering experiments 
were co-aligned mosaics of 3 samples to within half a 
degree with total masses of 3.21 g and 4.27 g for 3.5% 
and 7.5% Br samples respectively. Lattice parameters 
at 1.6 K were a = 7.86A, b = 6.96A, c = 6.08Aand an- 
gles a = 111.23°, /3 = 99.95°, 7 = 81.26° for 3.5% Br 
sample and a = 7.90A, b = 6.97A, c = 6.03A, and an- 
gles a = 111.15°, /3 = 99.91°, 7 = 81.29° for 7.5% Br 
sample. Time-of- flight neutron scattering experiments 
were carried out on the direct scattering geometry CNCS 
instruments^ at Spallation Neutron Source in Oak Ridge 
National Laboratory with 4.2 meV incident energy neu- 
trons. Triple axis neutron spectroscopy measurements 
were performed on TASP instrument at PSI.^° At TASP 
two modes of operation were used which we denote as 
the high resolution (hi-res) and the low resolution (lo-res) 
mode. In lo-res mode fixed 5 meV final energy and detec- 
tor distance of 75 cm with slit width of 30 mm was used. 
In hi-res mode fixed 3.5 me V final energy was used and 
detector analyzer distance was increased to 90 cm along 
with narrowing the slit to 20 mm in front of the detec- 
tor. Collimator of 80' was placed after the monochro- 
mator and liquid nitrogen cooled Be filter was employed 
between the sample and analyzer at all times. The reso- 
lutions determined from incoherent elastic scattering line 
widths at half maximum were 0.10 meV in hi-res mode 
and 0.25 meV in lo-res mode. Sample environment was a 
standard Orange cryostat in all experiments. 



III. RESULTS 

Low energy excitation spectrum up to 3meV was 
recorded by making full 360 degree rotations with 1° step 
at CNCS for 3.5% and 7.5% Br samples. Using scatter- 
ing plane determined by (1,0,0) and (0,1,1) reciprocal 
lattice vectors the 190° horizontal and 32° vertical ac- 
ceptance angles of the detector bank allowed probing of 
several Brillouin zones simultaneously, Fig{l] Recorded 



2 




-3-2-10123 
H, r.l.u. 

FIG. 1. (Color online) Top: Structure of PHCC, Magnetic 
Cu^^ (blue) ions are coupled via CI (green) ions to form 2D 
layers separated by piperazine molecules. Thick lines within 
layers indicate possible superexchange pathways. Bottom: 
Elastic neutron scattering intensity in the (H,L) plane (in- 
tegrated over K direction) for 3.5% Br sample showing the 
TOF experiment detector coverage. Thick red lines show lo- 
cations of triple axis scans projected onto (H,L) plane. 



events were projected to reciprocal sample coordinate 
system and binned to four-dimensional (H,K,L,E) matrix 
of size 101x101x101x96 spanning from (-3.5,-2.5,-2.5,- 
0.5) to (3.5,2.5,2.5,3.3) using Mantid^^ program. Two- 
dimensional projections from the full spectrum along the 
three main reciprocal directions are shown in Figures |2] 
and|3] In these figures the 6'(Q, uj) along the constant re- 
ciprocal coordinates was integrated over ±0.1 range and 
back-folded to first Brillouin zone. The dispersive gapped 
magnon mode is clearly visible in both samples. It is ev- 
ident that no drastic modification to the excitation spec- 
trum occurs in the bond disordered samples. Note that in 



the chosen sample orientation and instrument configura- 
tion we were unable to probe K = —L ^ and equivalent 
parts spectrum for E>0.3meV. 

Magnon linewidth and its wave vector dependence were 
studied in greater detail by triple axis spectroscopy in 
conjunction with well established resolution calculation 
methods using Popovici approximations^^ embedded in 
the Reslib package. Constant-Q scans at the zone center 
measured with hi-res mode of TASP at the temperature 
of 1.6 K were already shown in Fig. 4 of Ref.'S', As noted 
previously the intrinsic magnon linewidth and resonance 
energy are considerably increased with bond disorder. 
Although it is known from recent ESR studies^^ that 
the magnon band is split by 0.038 meV due to anisotropy, 
this splitting is small enough compared with our exper- 
imental resolution to justify using single mode approxi- 
mation (SMA). Morover, it has been established that for 
the 10%Br sample the splitting is reduced by a factor of 
two.'^ Hence the broadening is a lifetime effect induced 
by disorder. 

Temperature evolution of the magnon at zone center 
was studied on TASP using lo-res mode for the broader 
resonance of the 7.5% Br sample and in hi-res mode for 
the 3.5% Br sample. Constant-Q scans across magnon 
energy at different temperatures are shown in Fig. [6] 
Blue shift and damping with increasing temperature are 
evident for both samples. 

Triple axis spectroscopy was used to study the magnon 
wave vector dependence in lo-res mode at T=1.6K. Mea- 
sured spectra at different wave vectors were iteratively 
fitted with resolution convoluted Lorentzian lineshapes 
while adjusting the parameters in dispersion relation, 
FigsjS] and |4] Following the in-depth analysis of unper- 
turbed PHCCP we parametrized magnon dispersion: 

huj{Q) = (5o+7q)^ 

7Q = Bh cos(27r/i) + Bi cos(27r0 + Bhi[cos{27r{h + /)) + 
cos(27r(/i - /))] + B2h cos(47r/i) + B21 cos(47r/) (1) 

Best fitting parameters are collected in Table[l|for com- 
parsion. Fitted wave vector dependent magnon energies 
and linewidths are shown as symbols in Figj7| The key 
observations are that for Q=(0.5,0.5,-0.5) the gap in- 
creases from 1 meV to 1.36 and 1.50 meV for 3.5 and 
7. 5% Br samples respectively. At the same bandwidth is 
reduced by about 0.2 me V from 1.74 meV to 1.53 and 
1.58 meV. 

The bandwidth reduction is likely a consequence of 
modification of averaged exchange constants due of in- 
duced chemical pressure by Br substitution. In addi- 
tion we observed that magnon damping rate increases 
notably (by approximately a factor of three) when mov- 
ing towards zone boundaries. This drastic shortening of 
magnon lifetime will be discussed in the next section. 

Note that there exists small dispersion perpendicular 
to magnetic layers as seen from TOF spectra, middle 
panels of Figs. [2]and|3] However, in our fitting of disper- 
sion relation using triple axis spectroscopy data the used 
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FIG. 2. (color online) False color intensity maps of the dy- 
namic structure factor S{Q,(jj) for 3.5%Br PHCX sample 
along three main reciprocal lattice directions. In the above fig- 
ures the constant reciprocal coordinates denote center points 
of integration ranges over ±0.1 r.l.u. Data from distant unit 
cells has been back-folded to the first Brillouin zone. Note 
that some spurious scattering from instrumental origin is vis- 
ible at wave vectors close to direct beam. The solid lines are 
the fitted dispersion curves from triple axis data, Eq. [l]and 
Table. [l] The dashed line in the middle panel shows estimated 
interlayer bandwidth of 0.3 meV. 



sample orientation did not allow to extract the magni- 
tude of inter-layer coupling independently and hence we 
fixed it to zero as was done for unperturbed PHCC.'^ 
As the sample orientation was not optimized to measure 
dispersion along (0,1,0) direction a reliable fit cannot be 
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FIG. 3. (color online) False color intensity maps of the dy- 
namic structure factor S{Q,uj) for 7.5%Br PHCX sample 
along three main reciprocal lattice directions. Coordinates 
and lines as in Fig. |2] 



obtained from TOP spectra either. From the available 
data a crude estimate of ~0.3 meV for interlayer band- 
width can be given, shown by dashed lines in the middle 
panels of Figs. |2] and [3| The small interlayer coupling 
renormalizes the exact spin gap and bandwidth values 
slightly but the conclusions on the disorder effect remain 
qualitatively intact. 
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FIG. 4. (color online) Constant-Q scans across the magnon 
energy at different wave vectors (shown in parentheses) for 
the 7.5% Br sample. Corresponding solid lines are the 
fitted Lorentzian energy profiles convoluted with the four- 
dimensional instrumental resolution function. 



Parameter 


0% Br [13] 


3.5% Br 


7.5% Br 


in meV^ 








Bo 


5.44(2) 


5.50(4) 


6.12(7) 


B,. 


2.06(3) 


2.17(5) 


2.34(9) 


Bz 


1.07(3) 


1.02(5) 


1.28(8) 


B/iZ 


-0.39(1) 


-0.20(2) 


-0.05(4) 




-0.34(3) 


-0.15(4) 


-0.19(7) 


B2I 


-0.22(2) 


0.00(4) 


-0.04(6) 



TABLE I. Best fit parameters to describe the magnon dis- 
persion by Eq[l] for 3.5%Br and 7.5%Br PHCX samples in 
comparison with PHCC from Refll3l 



IV. DISCUSSION 

Random phase approximation (RPA) has been suc- 
cessful in describing the dispersion of crystal field ex- 
citations in several rare-earth compounds, as well as 
magnon dispersion in several dimerized transition metal 
based magnets. I^^M^ Even though the exact microscopic 
Hamiltonian remains unknown for PHCC and satisfac- 
tory coincidence with experimental spectrum has been 
obtained considering several magnetic Hamiltonian j^^ * ^^ 
the formal description of the spectrum Eq{l] is valid for 
RPA treatment. Namely, the gap is introduced by the 
effective dimer coupling Jq = \/Bo and the rest of the 
spectrum is dominated by smaller magnetic couplings 
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Hence we can write for the temperature re- 
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FIG. 5. (color online) Constant-Q scans across the magnon 
energy at different wave vectors (shown in parentheses) for 
the 3.5%Br sample. Corresponding solid lines are the fit- 
ted Lorentzian energy profiles convoluted with the four- 
dimensional instrumental resolution function. 



R{T) is the population difference of singlet ground state 
and the triply degenerate first excited state R{T) = 



(1 



-/3Jo^ 



)/(! + 36"^-^°). Intensity of the 



normalized gapP huj{T)Q = {Bq + R{T)jq)^ . Here 



ns - riT 

excitation is directly proportional to the population dif- 
ference the ground and excited state. Taking the inten- 
sity of the magnon peak at low temperature as scaling 
constant we expect I{T)/I{1.6K) = R{T). The finite 
temperature renormalization of magnon energy at dis- 
persion minimum the normalized intensity are shown in 
Figure [8) We see that for all samples the temperature 
induced blue shift and the intensity drop are explained 
by RPA treatment of an effective dimer model without 
any adjustable parameters. Temperature dependence of 
magnon linewidths will be discussed in detail elsewhere.!^ 
As a result of bond disorder the long lived magnon exci- 
tations get notably damped. Damping rate is smallest at 
the zone center Q=(0.5,0.5,-0.5) and increases by about a 
factor of three at zone boundaries for all measured direc- 
tions. This is in contrast to unperturbed PHCC where 
significant damping occurs at wave vectors where two- 



5 



Q = (0.5,0.5,-0.5) 
PHCX3.5%Br 



T-1.6K 





0.5 0.5 

(H 0.5 -0.5) (0.5 K-L) 



0.5 
(H 1+K- 



1-L) 



FIG. 7. (color online) Magnon excitation energies at 1.6 K 
and corresponding linewidths are shown with green circles 
and red squares for 3.5% and 7.5% Br samples respectively 
along three reciprocal space directions. Solid lines are the 
fitted dispersion curves described by Eq. [ijwith paremeters 
listed in Table |l] Shaded areas bounded by dotted, dashed 
and dash-dotted lines represent the two-magnon continua for 
7.5%, 3.5% and 0% Br samples, respectively. 



FIG. 6. (color online) Constant-Q scans across the magnon 
energy at different temperatures for the 3.5% and 7.5% Br 
samples. Corresponding solid lines are the fitted Lorentzian 
energy profiles convoluted with the four-dimensional instru- 
mental resolution function. 



magnon continuum has a minimum below single magnon 
spectrum.^ Similar albeit less rapid increase of damp- 
ing rate has been reported in Haldane spin chain com- 
pound CsNiClspSias opposed to IPA-CuCls where the sin- 
gle quasiparticle spectrum terminates altogether at the 
wave vector of the crossing.!^ Using the experimentally 
determined magnon dispersion relation, Eq{l] we calcu- 
late the energy of the lower boundary for the two-magnon 
continuum huj2L{Q) = min(ficj(Qi) + fia;(Q2))|Q2=Q-Qi 
and plot it as an envelopes of the shaded areas in Fig- 
ure [71 As a direct consequence of bandwidth reduction 
for PHCX samples the two magnon continuum drops be- 
low single particle spectrum only in a very narrow wave 
vector range, see the shaded areas in bottom right panel 
of Figl?! Hence, we can conclude that two-particle de- 
cay cannot account for the observed lifetime reduction in 
PHCX samples. 

On the other hand, in bond disordered samples the 
magnons can scatter off the impurities. Assuming this 
process is elastic, the available final state is a magnon 
with the same energy. In this case the scattering rate 
would follow the density of single magnon states (DOS) 



at particular energy T{uj) (x D{uj). To test this assump- 
tion we calculated D{(jj) numerically for both PHCX sam- 
ples by discretisizing the spectrum described by Eq{l]and 
convoluted it with the experimental energy resolution of 
the inelastic neutron experiment. These DOS curves are 
shown with solid line envelopes in Figj9j The comparison 
gives remarkably good qualitative agreement. This result 
highlights the difference between disorder and tempera- 
ture induced magnon lifetime shortening. In a disorder- 
free system magnon lifetimes shorten due to mutual colli- 
sions. However, in these collisional processes energy of a 
single magnon needs not be conserved and quali tatively 
different set of two magnon final states is allowed. I^SEUmI 



V. CONCLUSIONS 

We have measured the excitation spectrum up to 
4 meV for two bond disordered quasi-two-dimensional 
quantum magnets. No localized states were found to ap- 
pear within the spin gap and the excitation spectrum 
remains qualitatively similar to that of the unperturbed 
material. However, the spin gap increases and reduc- 
tion of bandwidth occurs. We attribute these effects to 
chemical pressure induced by Br substitution which af- 
fects the superexchange strengths. For both, 3.5 and 
7.5% Br samples RPA explains the temperature induced 
blue shift magnon mode at the zone centre. The drop 



6 




2 4 6 8 10 12 2 4 6 8 10 12 

T, K T, K 



FIG. 8. (Color online) Left: Gap renormalization as a func- 
tion of temperature: green circles and red squares are the 
magnon energies at Q=(0.5,0.5,-0.5) for 3.5% and 7.5%Br 
samples. Solid lines are the RPA predictions. Right: Inten- 
sity of the magnon peak normalized to 1.6 K intensity. Lines 
are population differences of singlet and triplet state, see dis- 
cussion in text. 
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FIG. 9. (color online) Intrinsic magnon linewidths for 3.5% 
Br (green circles) and 7.5% Br (red squares) PHCX samples 
collected from Fig]?! Dash-dot lines are the corresponding 
scaled densities of magnon states. Solid envelopes are the 
corresponding resolution convoluted densities of state. 

of magnon peak intensity with temperature is consistent 
with a singlet-triplet population difference of an effective 
dimer model. The long lived quasi-particles of the par- 
ent compound acquire finite lifetimes and the damping 
rate increases away from zone center in all directions. 
Damping rate in disordered samples follows the density 
of single magnon states. 
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